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Abstract 

MIMO systems are being considered as one of the key enabling technologies for future wireless 
networks. However, the decrease in capacity due to the presence of interferers in MIMO networks is 
not well understood. In this paper, we develop an analytical framework to characterize the capacity of 
MIMO communication systems in the presence of multiple MIMO co-channel interferers and noise. 
We consider the situation in which transmitters have no channel state information and all links undergo 
Rayleigh fading. We first generalize the determinant representation of hypergeometric functions with 
matrix arguments to the case when the argument matrices have eigenvalues of arbitrary multiplicity. This 
enables the derivation of the distribution of the eigenvalues of Gaussian quadratic forms and Wishart 
matrices with arbitrary correlation, with application to both single-user and multiuser MIMO systems. 
In particular, we derive the ergodic mutual information for MIMO systems in the presence of multiple 
MIMO interferers. Our analysis is valid for any number of interferers, each with arbitrary number of 
antennas having possibly unequal power levels. This framework, therefore, accommodates the study of 
distributed MIMO systems and accounts for different spatial positions of the MIMO interferers. 

Index Terms 

Eigenvalues distribution, Gaussian quadratic forms, Hypergeometric functions of matrix arguments, 
Interference, MIMO, Wishart matrices. 

I. Introduction 

The use of multiple transmitting and receiving antennas can provide high spectral efficiency 
and link reliability for point-to-point communication in fading environments [1], [2]. The analysis 
of capacity for MIMO channels in [3] suggested practical receiver structures to obtain such 
spectral efficiency. Since then, many studies have been devoted to the analysis of MIMO systems, 
starting from the ergodic [4] and outage [5] capacity for uncorrected fading to the case where 
correlation is present at one of the two sides (either at the transmitter or at the receiver) or at 
both sides [6]-[8]. The effect of time correlation is studied in [9]. 

Only a few papers, by using simulation or approximations, have studied the capacity of MIMO 
systems in the presence of cochannel interference. In particular, a simulation study is presented 
in [10] for cellular systems, assuming up to 3 transmit and 3 receive antennas. The simulations 
showed that cochannel interference can seriously degrade the overall capacity when MIMO links 
are used in cellular networks. In [11], [12] it is studied whether, in a MIMO multiuser scenario, 
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it is always convenient to use all transmitting antennas. It was found that for some values of SNR 
and SIR, allocating all power into a single transmitting antenna, rather than dividing the power 
equally among independent streams from the different antennas, would lead to a higher overall 
system mutual information. The study in [11], [12] adopts simulation to evaluate the capacity of 
MIMO systems in the presence of cochannel interference, and the difficulties in the evaluations 
limited the results to a scenario with two MIMO users employing at most two antenna elements. 
In [13] the replica method is used to obtain approximate moments of the capacity for MIMO 
systems with large number of antenna elements including the presence of interference. The 
approximation requires iterative numerical methods to solve a system of non-linear equations, 
and its accuracy has to be verified by computer simulations. A multiuser MIMO system with 
specific receiver structures is analyzed for the interference-limited case in [14], [15]. 

The MIMO capacity at high and low SNR for interference-limited scenarios is addressed in 
[16], [17]. A worst-case analysis for MIMO capacity with CSI at the transmitter and at the 
receiver, conditioned on the channel matrix, can be found in [18]. Asymptotic results for the 
Rician channel in the presence of interference can be found in [19]. 

In this paper, we develop an analytical framework to analyze the ergodic capacity of MIMO 
systems in the presence of multiple MIMO cochannel interferers and AWGN. Throughout the 
paper we consider rich scattering environments in which transmitters have no CSI, the receiver 
has perfect CSI, and all links undergo frequency flat Rayleigh fading. The key contributions of 
the paper are as follows: 

• Generalization of the determinant representation of hypergeometric functions with matrix 
arguments to the case where matrices in the arguments have eigenvalues with arbitrary 
multiplicity. 

• Derivation, using the generalized representation, of the joint p.d.f. of the eigenvalues of 
complex Gaussian quadratic forms and Wishart matrices, with arbitrary multiplicities for 
the eigenvalues of the associated covariance matrix. 

• Derivation of the ergodic capacity of single-user MIMO systems that accounts for arbitrary 
power levels and arbitrary correlation across the transmitting antenna elements, or arbitrary 
correlation at the receiver side. 

• Derivation of capacity expressions for MIMO systems in the presence of multiple MIMO 
interferers, valid for any number of interferers, each with arbitrary number of antennas 
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having possibly unequal power levels. 
The paper is organized as follows: in Section II we introduce the system model for multiuser 
MIMO setting, relating the ergodic capacity of MIMO systems in the presence of multiple 
MIMO interferers to that of single-user MIMO systems with no interference. General results 
on hypergeometric functions of matrix arguments are given in Section III. The joint p.d.f. of 
eigenvalues for Gaussian quadratic forms and Wishart matrices with arbitrary correlation is given 
in Section IV. In Section V we give a unified expression for the capacity of single-user MIMO 
systems that accounts for arbitrary correlation matrix at one side. Numerical results for MIMO 
relay networks and multiuser MIMO are presented in Section VI, and conclusions are given in 
Section VII. 

Throughout the paper vectors and matrices are indicated by bold, |A| and det A denote the 
determinant of matrix A, and a it j is the element of A. Expectation operator is denoted 

by E{-}, and in particular Ex{-} denotes expectation with respect to the random variable X. 
The superscript f denotes conjugation and transposition, I is the identity matrix (in particular I n 
refers to the (n x n) identity matrix), tr {A} is the trace of A and © is used for the direct sum 
of matrices defined as A © B = diag (A, B) [20]. 



We consider a network scenario as shown in Fig. 1, where a MIMO-(A r T0 , N R ) link, with iV T0 
and N R denoting the numbers of transmitting and receiving antennas, respectively, is subject to 
Ni MIMO co-channel interferers from other links, each with arbitrary number of antennas. The 
V R -dimensional equivalent lowpass signal y, after matched filtering and sampling, at the output 
of the receiving antennas can be written as 



where x , xi, . . . , xjy, denote the complex transmitted vectors with dimensions Nj , N T1 , . . . , N TNl , 
respectively. Subscript is used for the desired signal, while subscripts 1, . . . , Nj are for the 
interferers. The additive noise n is an V R -dimensional random vector with zero-mean i.i.d. 
circularly symmetric complex Gaussian entries, each with independent real and imaginary parts 



II. System Models 




(1) 



k=l 



having variance a /2, so that Ejnn 1 "} 



<j 2 I. The power transmitted from the k- 



•— user is 




DRAFT 



4 



SUBM. TO IEEE TRANS. ON INF. TH. 



Po, x 



N- 



< 
< 



TO 



Pi, X X 



Pn v x JVi 



< 



TJVi 



H r 



A 



H 



11 



Fig. 1. MIMO Network. 



The matrices H fc in (1) denote the channel matrices of size (iV R x N Tk ) with complex elements 
hfj describing the gain of the radio channel between the j— transmitting antenna of the k— 
MIMO interferers and the i— receiving antenna of the desired link. In particular, H is the 
matrix describing the channel of the desired link (see Fig. 1). 

When considering statistical variations of the channel, the channel gains must be described as 
r.v.. In particular, we assume uncorrelated MIMO Rayleigh fading channels for which the entries 
of Hfc are i.i.d. circularly symmetric complex Gaussian r.v. with zero-mean and variance one, i.e., 
E jj/i^j | 2 | = 1. With this normalization, Pk represents the short-term average received power 
per antenna element from user k, which depends on the transmit power, path-loss, and shadowing 
between transmitter k and the (interfered) receiver. Thus, the Pk are in general different. 

Conditioned to the channel matrices {H fc }^ , the mutual information between the received 
vector, y, and the desired transmitted vector, x , is [21] 

X (x ; y | {H fc }£ ) = H (y | {H fc }fi ) - H (y | x , {H fe }£ ) (2) 

DRAFT 



CHIANI, WIN, SHIN: MIMO NETWORKS. 



5 



where TC(-) denotes differential entropy. 

Here we consider the scenario in which the receiver has perfect CSI, and all the transmitters 
have no CSI. Note that the term CSI includes the information about the channels associated with 
all other MIMO interfering users. In this case, since the users do not know what is the interference 
seen at the receiver (if any), a reasonable strategy is that each user transmits circularly symmetric 
Gaussian vector signals with zero mean and i.i.d. elements. Thus, the transmit power per antenna 
element of the k— user is Pk/N Tk . Note that this model includes the case in which the power 
levels of the individual antennas are different: it suffices to decompose a transmitter into virtual 
sub-transmitters, each with the proper power level. 

Hence, conditioned on all channel matrices {H fc }^ in (1), both y and y|x are circularly 
symmetric Gaussian. Since the differential entropy of a Gaussian vector is proportional to 
the logarithm of the determinant of its covariance matrix, we obtain the conditional mutual 
information 



where K y and K y | X(J respectively denote the covariance matrices of y and y|x , conditioned 
on the channel gains {Hfc}^ . By expanding the covariance matrices using (1), the conditional 
mutual information of a MIMO link in the presence of multiple MIMO interferers with CSI 
only at the receiver is then given by: 




(3) 




(4) 



where the N R x (X)*=i N Ti ) matrix H is 



H = [Hi|H 2 | • ■ ■ [HjvJ 



the N R x (Y,i=o N ?i) matrix H is 



H= [H |H] 



the covariance matrices \I> are 



* = Qt 1 Nji © g 2 In T2 © ■ • • © QNi IjVi 



(5) 



and 



* = ft) IiV T0 © * 



(6) 
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with 

Pi 

61 ~ N Ti a 



(7) 



With random channel matrices the mutual information in (4) is the difference between random 
variables of the form logdet (I + H3>H^) where the elements of H are i.i.d. complex Gaussian 
and $ is a covariance matrix. The statistics of such random variables has been investigated 
in [6]-[8], assuming that the eigenvalues of <fr were distinct. However, in the scenario under 
analysis these results cannot be used directly, since in (4) each eigenvalue Qi of ^ and *Sf has 
multiplicity N Ti . 

We consider the ergodic mutual information as a performance measure: taking the expectation 
of (4) with respect to the distribution of {H fc }^ , we get 

C MU = E {CW ({H fe }fi ) } = C su N Tl , N R , * J - C su N Tl , N R , * j (8) 

where Csu (n-T,n R ,<f>) — Eh{ logdet (l„ R + H$H^)} denotes the ergodic mutual information 
of a single-user MIMO-(n T , n R ) Rayleigh fading channel with unit noise variance per receiving 
antenna and channel covariance matrix $ at the transmitter. 

Note that the "building block" E H {log det (I + H$H^) } is simple to evaluate when the 
covariance matrix <fr is proportional to an identity matrix, which corresponds to a typical 
interference-free case with equal transmit power among all transmitting antennas (see, e.g., 
[4]). In contrast, in the presence of interference, the covariance matrix is of the type indicated 
in (5) and (6), where the power levels of the different users are in general different. Note that 
even when the power for the i— user is equally spread over the Nj { antennas, the matrices 
in (5) and (6) are in general not proportional to identity matrices and their eigenvalues have 
multiplicities greater than one. Therefore, studying MIMO systems in the presence of multiple 
MIMO cochannel interferers requires the characterization of Csu {nr, n R , <&) in a general setting 
in which the covariance matrix $ has eigenvalues of arbitrary multiplicities. 

To this aim, we derive in the next sections simple expressions for the hypergeometric functions 
of matrix arguments with not necessarily distinct eigenvalues; then, we obtain the joint p.d.f. 
of the eigenvalues of central Wishart matrices as well as that of Gaussian quadratic forms with 
arbitrary covariance matrix. 
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III. Hypergeometric functions with matrix arguments having arbitrary 

EIGENVALUES 



Hypergeometric functions with matrix arguments [22] have been used extensively in multi- 
variate statistical analysis, especially in problems related to the distribution of random matrices 
[23]. These functions are defined in terms of a series of zonal polynomials, and, as such, they 
are functions only of the eigenvalues (or latent roots) of the argument matrices [22], [23]. 

Definition 1: The hypergeometric functions of two Hermitian m x m matrices A and W are 
defined by [22] 



,F q (ai,...,a p , b q] A, W) = g £ ^ . . . ^ (9) 

where C K (-) is a symmetric homogeneous polynomial of degree k in the eigenvalues of its 
argument, called zonal polynomial, the sum ^ K is over all partitions of k, i.e., k — (ki, . . . , k m ) 
with k\ > k 2 > ■ ■ ■ > k m > 0, k\ + ki + ■ • • + k m = k, and the generalized hypergeometric 
coefficient (a) K is given by (a) K = YliLi { a ~ — 1))^. w ^ m ( a )fc = a ( a + 1) • • • (a + A; — 1), 
(a) = 1. 

We remark that zonal polynomials are symmetric polynomials in the eigenvalues of the matrix 
argument. Therefore, hypergeometric functions are only functions of the eigenvalues of their 
matrix arguments. In other words, without loss of generality we can replace A and W with the 
diagonal matrices diag (Ai, . . . A m ) and diag (wi, . . . w m ), where Aj and Wj are the eigenvalues 
of A and W, respectively. Clearly the order of A and W is unimportant. 

It is quite evident that these functions expressed as a series of zonal polynomials are in general 
very difficult to manage and the form of (9) is not tractable for further analysis. Fortunately, 
when the eigenvalues of A and W are all distinct, a simpler expression in terms of determinants 
of matrices whose elements are hypergeometric functions of scalar arguments can be obtained 
as follows [24, Lemma 3]: 

Lemma 1: ([Khatri, 1970]) Let A = diag (Ai, . . . A m ) and W = diag (wi, . . . w m ) with Ai > 
• • • > A m and w\ > ■ ■ ■ > w m . Then we have 

^ (m \b) IGI 
t F, (oi , .... a r , K ... , K A, W) = ^(m) 4m){a)nj<j{K _ Xj)Ui<Awi _ aj) (10) 
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where F 



(m) 



n 



i + I) 1 1 and the ij— element of 



the (to x to) matrix G is defined in terms of hypergeometric functions of scalar arguments as 
follows 

,a p ;6i,...,6 g ;AiWj) (11) 



9u 



pFq [ a l, 



with hi = di — m + 1 and bi — bi — to + 1. 
Important particular cases are 



F (A,W) =r 



IGr 



(m) 



TO 



and 



iF (r;A,W) 



(m) 



rn 



Ui<j (Ai - Aj) (^i - ^i) 
I Gil 



^! m) (r) n,<i - X j) Ui<j (*>i - v>i) 



(12) 



(13) 



where the ij— elements of Go and Gi are given by e XiWj and (1 — A 



* W 3 



i m—r—l 



, respectively. 



These expressions have been recently used to study the distribution of Gaussian quadratic 
forms, to express the p.d.f. of the eigenvalues of Wishart matrices, and to analyze the information- 
theoretic capacity and error rates of communication systems involving multiple antennas [5]-[8], 
[25]— [31]. However, it is important to underline that Lemma 1 requires the eigenvalues of the 
matrices to be all distinct. 

Here, we generalize Lemma 1 to include the case where the eigenvalues are not necessarily 
distinct. To this aim we first need the following lemma. 

Lemma 2: Let P : A — > R be defined over A C R m as follows: 



P(wi, ...,w r . 



n 



, l<j {w l -w j/ 



AM AM ■•■ AM) 



(14) 



fm{Wl) f m {w 2 ) ■ ■ ■ fm( 

where w\ > w 2 ■ ■ ■ > w m , and the functions fi(w) have derivatives f- (w) = d of orders 
at least to — 1 throughout neighborhoods of the points w 1 , . . . , w m . 

Then, the continuous extension P(wi, io 2 , • • • , w m ) of the function P(wi,W2, ■ ■ ■ ,w m ) to 
those points in R m with L coincident arguments wk = wk+i = is obtained by 
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removing the zero factors from the denominator in (14), replacing the columns of the matrix 
in (14) corresponding to the coincident arguments with the successive derivatives fl L ~ l \wK), 
I = 1, . . . , L, and dividing by a scaling factor F(l)(L) = Yli=i 
For example, for W\ = w 2 = ■ ■ ■ wl, this procedure gives 

1 



P(w 1 ,w 2 , 



W r , 



f[ L - l \ Wl ) ff- 2 \ Wl ) ■ ■ ■ f x { Wl ) h{w L+x ) ■■■ h{w m ) 



(15) 



fm 1} Ol) frt 2) Ol) ••• f m {Wl) f m {w L +l) fm{w m ) 

More generally, a similar expression is valid if there are more groups of coinciding arguments: 
in this case, for each group of coincident arguments wk = . . . = wk+l-i the correspondent 
columns of the matrix in (14) are to be replaced by f- L ~ l \wK), I = 1,...,L, with a scaling 
factor Yli=i 



Proof: See Appendix I. 

With Lemma 2 we can now generalize (12), (13) and, more generally, (10). 
Lemma 3: Let A = diag (Ai, . . . A m ) and W = diag (wi, . . . w m ) with Ai > 

Wk+L-i > W k+L > ■ ■ ■ > w m . Then we have 1 
IV) (m) |G| 



□ 



> A m and 



w x > ■ ■ ■ > W k = W k+ i = ■ ■ 
F (A, W 



r (L){L) Ui<j (Ai - Xj) II,. K - Wj) 



(16) 



where the elements of G are 



X L-l+k-j eXi 



9i,j 



(17) 



j — k, . . . , k + L — 1 
e XlWj elsewhere 
that is, the matrix G is the same as that appearing in (12) except that the L columns corresponding 
to the coincident eigenvalues are Af - 1 e XiWk , Xf ~ 2 e x ' Wk , . . . , Xfe x ' Wk , Xie x > Wk , e XlWk . 

Proof: The proof is immediate by direct application of Lemma 2 with fiiyS) = e XlW . 

□ 



From here on we will use the same symbols for the functions (12), (13), (10) and their continuous extension. 
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Lemma 3 can be directly extended to more groups of coincident eigenvalues. In general, 
the rule is that each eigenvalue w of multiplicity L > 1 gives rise to L columns Af _1 e Ai,u , 
\f~ 2 e XiW , . . . , Xje XlW , A;e Al ™, e XiW in the matrix G of (16), with the proper scaling factor r (L) (L). 

Using Lemma 3 with k = m — L + 1 and Wk = results in the following corollary, for the 
case where some eigenvalues are equal to zero. 

Corollary 1: Let A = diag (Ai, . . . X m ) and W = diag (wi, . . . w m ) with Ai > • • • > A m and 
Wi> • • - > w m ^L+i = w m ^L + 2 = ■ ■ ■ = w m = 0. Then we have 

oF (A, W) = 54^1 B (18) 

where the elements of G are as follows 

J K~ 3 j = m- L + l,...,m 
9i,j = < . (19) 
I e lU "-> elsewhere. 

We can apply a similar methodology to derive the general expression for iFq(-; •, •), as in the 
following Lemma. 

Lemma 4: Let A = diag (Ai, . . . A m ) and W = diag (wi, . . . w m ) with Ai > • • • > A m and 
W! > ■ ■ ■ > w k = w k+ i = ■ ■■ = Wk+L-i > w k+ L > ■■■ > w m . Then we have 

F (v AW)- r HW(-l) W 7 L ~ 1 (7-l) L - 2 ---(7-^ + 2) |A| 

where 7 = m — r — 1 and the (m x m) matrix A has elements as follows 

/ Af- W (1 - X lWj r iL - 1+k ^ j = k, . . . , k + L - 1 
= < (21) 
I (1 — XiWj) J elsewhere. 

In other words, the matrix A is the same as that appearing in (13), except that the L columns cor- 
responding to the L coincident eigenvalues are Af _1 (1 — Ajitffc)" 1 '~^ L ~ 1 > ; . . . , Aj (1 — AjU^) 7-1 , 
(1 - XiWkf. 

Proof: For the proof we apply Lemma 2 with fi(w) = (1 — XiWk) 1 , whose n— derivative 

is fl n \w) = (-A,)" 7 (7 - 1) • " " (7 - n + 1) (1 - X lWk y- n . □ 

Lemma 4 can be further generalized to more groups of coincident eigenvalues: each eigenvalue 
w of multiplicity L > 1 gives rise to L columns Xf^^-Xiw)^ 1 '^, . . . , Af (l-Ajw) 7-2 , A;(l- 
Aitu) 7-1 , (1 - Xiw)~t in the matrix A of (20), and to a factor (-l)( L - 1 ) L / 2 ~ f L ~ 1 ■ • • (7 — L + 
2)/T ( L)(i). 
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Using Lemma 4 with k = m — L + 1 and w k = results in the following corollary. 
Corollary 2: Let A = diag (Aj., . . . A m ) and W = diag (wx, . . . w m ) with Ai > • • • > A m and 
»!>■••> w m _ L+l = w m _ L+2 = ■ ■ ■ = w m = 0. Then we have that (20) holds, with 

. K~ 3 j = in — L + 1, . . . ,m 

a i,j = \ . ( 22 ) 



, 1 — XiWj) 1 elsewhere. 



In other words, the matrix A has in this case the last L columns with elements Xf ,Xf , . . . , A i; 1. 
Finally, we give the result for the p F q (-). 

Lemma 5: Let A = diag (Ai, . . . X m ) and W = diag (wi, . . . w m ) with Ai > • • • > A m and 

Wx > ■ ■ ■ > w h = w k+1 = ■■■ = w k+L ^i > w k+L > ■■■ > w m . Then we have 

|C| 

p F q (ax, . . . ,a p ;bx,. . . , b q ; A, W) = S ^ r — r (23) 

where the (m x m) matrix C has elements as follows 



Xf 1+k j p F q (ax - m + L + k - j, . . . ,b q - m + L + k - j; XiWj) (24) 
for j = k, . . . , k + L — 1, and 

elsewhere. In (23) the constant S is 

_ = r (m) (m) ipj m) (b) yr (oi)i(o 2 )i • • • (a P )i 
a ~ T (L) (L) 4 m \a)l\ (bxHb2)i---(b q )i' 
Proof: See Appendix I. □ 

IV. Gaussian quadratic forms with covariance matrix having eigenvalues of 

ARBITRARY MULTIPLICITY 

We now derive the joint p.d.f. of the eigenvalues for Gaussian quadratic forms and central 
Wishart matrices with arbitrary one-sided correlation matrix. 

Lemma 6: Let H be a complex Gaussian (pxn) random matrix with zero-mean, unit variance, 
i.i.d. entries and let $ be an (n x n) positive definite matrix. The joint p.d.f. of the (real) non-zero 
ordered eigenvalues Ai > A 2 > . . . > A nmin > of the (p x p) quadratic form W = H$H^ is 
given by 

^min 

f x (xx,...,x nn J = K\V(x)\ G(x,/x) n*r nmm (25) 



i=l 
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where n m j n = min(ra, p), V(x) is the (n m i n x n m [ n ) Vandermonde matrix with elements Vij = x* , 



(_l)p(n-n rain ) nf=l^f 



r (n m ,„)(p) n£=i r (tn i )(m i )ni < ,-(^«)-^(j)) miny ( } 

and > /j(2) . . . > /i(L) are the L distinct eigenvalues of <fr _1 , with corresponding multiplicities 
mi, ... , mi such that J2i=i m « = n - 

The (n x n) matrix G(x, fi) has elements 

- /(-ar^e-^ j = l,...,n mm 

9i,j - S __ d (27) 

[ [n-fidiVieJ 1 J =n min + l,...,n 

where [a] k = a(a — 1) ■ ■ ■ (a — k + 1), [a] = 1, denotes the unique integer such that 

mi + . . . + m e ._i < i < rrii + . . . + m ei 

and 



ei 

di = }] m fc - i. 



k=i 

Proof: See Appendix I. □ 
Note that Lemma 6 gives, in a compact form, the general joint distribution for the eigenvalues 

of a central Wishart (p > n), and central pseudo-Wishart or quadratic form (n > p), with 

arbitrary one-sided correlation matrix with not-necessarily distinct eigenvalues. 

In fact, Lemma 6 can be used for both p > n and n > p; in particular, for n > p we have 

nr=i x "i~ nmm = 1 in (25), while for p > n the second row in (27) disappears and (— l)p( n_n mm) — ]_ 

in (26). 

Moreover, using Lemma 6 and the results in [32], [33] we can also derive the marginal 
distribution of individual eigenvalues or of an arbitrary subset of the eigenvalues. 

V. Ergodic Mutual information of a single-user MIMO system 

In this section we provide a unified analysis of the ergodic mutual information of a single-user 
MIMO system with arbitrary power levels/correlation among the transmitting antenna elements 
or arbitrary correlation at the receiver, admitting correlation matrices with not-necessarily distinct 
eigenvalues. 

Let us consider the function 

C su (n,p, $) = E H {logdet (l p + H^H 1 ") } (28) 
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where $ is a generic (n x n) positive definite matrix and H is a (p x n) random matrix with 
zero-mean, unit variance complex Gaussian i.i.d. entries. 

Now, consider a single-user MIMO-(ra T , n R ) Rayleigh fading channel with \E , t, ^ r denoting 
the (n T x rij) transmit and (n R x n R ) receive correlation matrices, respectively, having diagonal 
elements equal to one. Assume the transmit vector x is zero-mean complex Gaussian, with 
arbitrary (but fixed) (n T x n T ) covariance matrix Q = Ejxx^} so that tr{Q} = P. Then, 
the function (28) can be used to express the ergodic mutual information in the following cases 



1) the MIMO-(n T , n R ) channel with no correlation at the receiver (^ R = I), covariance matrix 
at the transmitter side ^t, and transmit covariance matrix Q. 

In this case the mutual information is Csu (n T ,n R ,Q>) with = (1/°" 2 )*tQ- If also 
* T = I, we have <E» = (l/a 2 )Q and therefore tr{$} = P/a 2 . 

2) the MIMO-(n T , n R ) channel with no correlation at the transmitter (^ T = I), covariance 
matrix at the receiver side \I/ R , and equal power allocation Q = P/n T l. 

In this case the capacity is Csv (n R ,n T ,^) with ^ = (-PMt0" 2 ) , I , r> giving tr{$} = 
(P/a 2 )(n R /riT), in accordance to [6, Theorem 1]. 

In both cases P/a 2 represents the SNR per receiving antenna. 

By indicating with n m { n = mm(n,p) and with /a(-, •••,■) the joint p.d.f. of the (real) ordered 
non-zero eigenvalues Ai > A2 > . . . > A„ min > of the (p x p) random matrix W = H$H^, 
we can write: 



where the multiple integral is over the domain V md = {00 > x\ > x 2 > . . . > x„ mm > 0} and 

dx = dxi dxo ■ ■ ■ dx n . . 

The nested integral in (29) can be evaluated using the results from previous sections and 
Appendix II, leading to the following Theorem. 

Theorem 1: The ergodic mutual information of a MEMO Rayleigh fading channel with CSI at 
the receiver only and one-sided correlation matrix $ having eigenvalues of arbitrary multiplicities 



[6]-[8]: 





i=l 



(29) 
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is given by 

C su (n,p,$)=K^det(R( fe )). (30) 

k=l 

In the previous equation n min = min(n,p), the matrix has elements 



io°° x p ~ nmm+j ~ l+dl e- x ^dx j = 1, . . . , n min , j ^ Jfe 

Jo°° sP-^^'-^e- 1 log (l + x)dx j = 1, . . . , n min , j = jfe (31) 

k [n - ik J = rimm + 1, . . . , n 

and [a]fc, e*, dj, K are defined as in Lemma 6, where fj,m > fi/2) ■ ■ ■ > Ml) are the L distinct 
eigenvalues of with corresponding multiplicities mi, . . . , m^. 

In Section IV it is shown that the joint p.d.f. of the ordered eigenvalues of W can be written 
as (25), where the elements of V(x), G(x, //) are real functions of xi, . . . , x nmin . Thus, by using 
Appendix II, the multiple integral in (29) reduces to (30). □ 

Note that the integral in (31) can be evaluated easily with standard numerical techniques; 
however, the integral can be further simplified, using the identities x m e~ xll dx = m!//i m+1 , 
and J °° x m e~ x ^ ln(l + x)dx = m!e M ^™ o r(i — m, where •) is the incomplete 

Gamma function. 

Theorem 1 gives, in a unified way, the exact mutual information for MIMO systems, encom- 
passing the cases of n R > n T and n T > n R with arbitrary correlation at the transmitter or at 
the receiver, avoiding the need for Monte Carlo evaluation. The application of the results in 
Sections III-V enables a unified analysis for MIMO systems, which allow the generalization 
for ergodic and outage capacity [6]-[8], [29], for optimum combining multiple antenna systems 
[26], [27], for MIMO-MMSE systems [28], for MIMO relay networks [34], [35], as well as for 
multiuser MIMO systems and for distributed MIMO systems, accounting arbitrary covariance 
matrices. For example, after the first derivation of the hypergeometric functions of matrices with 
non-distinct eigenvalues in [36], other applications to multiple antenna systems have appeared 
in [32], [37]-[40]. 



VI. Numerical Results 
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Let us first apply Theorem 1 to the analysis of a single-user MEMO system with unequal 
power levels among the transmitting antennas. Figure 2 shows the ergodic mutual information 2 
of a MIMO-(6, 3) Rayleigh channel, where the relative transmitted power levels are {1 + A, 1 + 
A, 1 + A, 1 — A, 1 — A, 1 — A}. The particular cases A = and A = 1 are equivalent to the 
equal power levels over 6 and 3 transmitting antennas, respectively. This figure shows how the 
capacity decreases as A increases from to 1, with a behavior in accordance to analysis based 
on majorization theory [41]. 

As another example of application, we evaluate the performance of MEMO relay networks in 
Rayleigh fading [34], [35]. For such networks the network capacity is upper bounded by [35, 
eq. (5)], which can be easily put in the form C u = lEnjlogdet (I + H<3?H^)}, and evaluated 
in closed form by Theorem 1. In Fig. 3 we report the exact C u as obtained from Theorem 1, 
compared with the Jensen's inequality [35, Theorem 1]. The figure has been obtained for a 
source node with 4 antennas, 5 relays each equipped with 2 antennas, as a function of the total 
equivalent SNR here defined as SNR = tr {<!>}. We assume, for the 5 relays, that the received 
power is distributed proportionally to the weights {1, 2, 5, 10, 20}. It can observed that the results 
based on the Jensen's inequality can be overly optimistic. 

As a third example of application we evaluate, using (8) together with Theorem 1, the exact 
expression of the ergodic mutual information of MIMO systems in the presence of multiple 
MIMO interferers in Rayleigh fading. In particular, the eigenvalues to be used in Theorem 1 are 
given by ixia = 1/ Qi = a 2 N Ti /Pi, allowing an easy analysis for several scenarios. We define 
the average SNR per receiving antenna as SNR = Po/a 2 giving g = SNR/iV T0 , and the SIR 
as SIR = PqJ Ylii>\ Pi- 3 Fig- ^ shows the ergodic mutual information for a MIMO-(6, 6) system 
as a function of the SIR, in the presence of one MIMO cochannel interferer having iV T1 equal 
power transmitting antennas. It can be noted that the capacity decreases with the increase in the 
number of interfering antenna elements, tending to the curve obtained by using the Gaussian 
approximation. 4 Despite the fact that the received vector y in (1) is Gaussian conditioned on 

2 For the numerical results we use the base 2 of logarithm in all formulas, giving a mutual information in bits/s/Hz. 

3 We recall that, with our normalization on the channel gains, the mean received power from user i is P i; and our definition 
of SIR account for the total interference power. 

4 With Gaussian approximation the performance is evaluated as if interference were absent, except the overall noise power is 
set to a 2 + X^i>i Pii g' vm 8 a signal-to-interference-plus-noise ratio SINR = + gjp) 1 . 
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the channel matrices, and that the elements of are Gaussian, approximating the cumulative 
interference as a spatially white complex Gaussian vector is pessimistic for analyzing MIMO 
systems in the presence of interference, unless the number of transmitting antenna of the interferer 
is large compared with that of the desired user. This is because the Gaussian approximation 
implicitly assumes that the receiver does not exploit the CSI of the interferers (single-user 
receiver), whereas the exact capacity accounts for the knowledge of all CSI at the receiver. In 
the same figure we also report, using circles, the capacity of a single-user MIMO-(iV T0 , N R —N T1 ) 
for N R > iV T1 . It can be observed that the capacity of the MIMO-(A^ T0 , N R ) in the presence of 
Nji interfering antenna elements approaches asymptotically, for large interference power, to a 
floor given by the capacity of a single-user MIMO-(iV T0 , N R — N T1 ) system. This behavior can 
be thought of as using Nj 1 DoF at the receiver to null the interference in a small SIR regime. 
On the other hand, when N R < N T1 the capacity approaches to zero for small SIR. This is due 
to the limited DoF at the receiver (related to the number N R of receiving antenna elements) that 
prevents mitigating all interfering signals (one from each antenna elements) while, at the same 
time, processing the iV T0 useful parallel streams, as previously observed for multiple antenna 
systems with optimum combining [2], [26], [27]. 

Finally, in Fig. 5 we consider a MIMO-(A r T0 , 6) system in the presence of one and two 
MIMO interferers in the network, each equipped with the same number of antennas as for the 
desired user. We clearly see here two different regions: for small SIR the interference effect is 
dominant, and it is better for all users to employ the minimum number of transmitting antennas 
(i.e., MIMO-(3, 6) for all users), so as to allow the receiver to mitigate the interfering signals. 
On the contrary, for large SIR the channel tends to that of a single-user MIMO system and it 
is better to employ the maximum number of transmitting antennas. In the same figure we also 
report the capacity for interference-free channels, which represents the asymptotes of the four 
curves, as well as the Gaussian approximation, which incorrectly indicates that it is always better 
to use the largest possible number of transmitting antennas. 

It can be also verified that, in a network where all nodes are using the same MIMO-(n, n) 
systems, larger values of n achieve higher mutual information, for all values of SIR and SNR. 
Note, however, that when increasing the number of antennas and users, correlation may arise in 
the channel matrices. 
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VII. Conclusion 

We have studied MIMO communication systems in the presence of multiple MIMO interferers 
and noise. To this aim, we first generalized the determinant representations for hypergeometric 
functions with matrix arguments to the case where the eigenvalues of the argument matrices have 
arbitrary multiplicities. Then, we derived a unified formula for the joint p.d.f. of the eigenvalues 
for central Wishart matrices and Gaussian quadratic forms, allowing arbitrary multiplicities for 
the covariance matrix eigenvalues. These new results enable the analysis of many scenarios 
involving MEMO systems. For example, we derived a unified expression for the ergodic mutual 
information of MIMO Rayleigh fading channels, which applies to transmit or receive correlation 
matrices with eigenvalues of arbitrary multiplicities. We have shown how to apply the new 
expressions to MIMO networks, deriving in closed form the ergodic mutual information of 
MIMO systems in the presence of multiple MIMO interferers. 



Appendix I 
Proofs 

A. Proof of Lemma 2 

For ease of notation and without loss of generality we consider the case K = 1, where the 
application of the lemma leads to (15). For the proof we proceed by induction. First, the result 
in (15) is obvious for L — 1, since in this case (15) coincides with (14). Then, we must show 
that if (15) is true for any L then it is also true for L + 1. So, assuming that (15) holds for L, 
we must find 

lim P(w 1 ,...,w m ). 

In this regard note that, with w% = w 2 = ■ ■ ■ = w L the product Yl%<j w^w ■ ( w i ~ w j) m (1^) 
contains exactly L factors with value e = wl — Wl+i- Thus, by rewriting wl+i = wl — e we 
have 



lim P(w 1 ,...,w m ) = — 

UiKj^^L+M ~ Wj) Ui=l 
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■ lim — ; 



f[ L 1] {wl) ■■■ fi{w L ) fi(w L -e) ■■■ fi(w m ) 



fm 1} (W L ) ■ ■ ■ f m {w L ) f m (w L - e) • • • f m (w m ) 



We can now apply the Taylor expansion to the functions 

L 



w 



n=0 



-e 



nl 



+ 0(e 



L+U 



(32) 



(33) 



where 0(e) denotes omitted terms of order e. We also know from basic algebra that, seen as a 
function of a column with the others fixed, the determinant is a linear function of the entries in 
the given column, as is clear for example from the Laplace expansion. Therefore, we have 



lim P(wi, . . . , w Ti 



n 



i<j,Wi^Wj,j^L+l 



Wi 



-e 



1 L 

e L ^— ' n 

n=0 



/i (L_1) K) ■•• AK) /! n) K) ■■■ AW 



f!n ^{wl) ■■■ fm{w L ) fm'{w L ) ••• f m (w m ) 



in). 



(34) 



In the summation above the determinants for n = 0, . . . , L — 1 are zero since there are 
coincident columns. So, in the limit for e — > only the term of grade L remains. 

By simplifying and reordering the first L + 1 columns of the matrix in (34), with a cyclic 
permutation having sign equal to (— 1) L , we finally have 



lim P(wi,...,w % 



1 



Y[i<j t v,^wX w i- w i)Xii=\ 



'/.! 
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fi K+i) ••• /iK + i) h(w L+2 ) ■■■ fi(w m ) 



(35) 



fm'(w L+1 ) ■■■ f m (w L+1 ) f m {w L+2 ) ■■■ f m (w m ) 

which is again in the form of (15). This concludes the proof by induction of Lemma 2 for 

w\ = ■ ■ ■ = w L . 

The extension to different K and more groups of coincident arguments is straightforward. 

B. Proof of Lemma 5. 

The derivatives of the hypergeometric function of scalar arguments can be expressed as 

— pF q (a 1 , . . . , a p ; 6i, . . . , b q ; z) = — 7r^ P F q (°i + n , ■ ■ ■ , a v + "5 h + n, . . . , b q + n; z) . 



gives Lemma 5. 

C. Proof of Lemma 6 

Here, based on Section III, we prove Lemma 6 concerning the eigenvalues distribution of 
Gaussian quadratic forms. The problem is related to the distribution of random matrices of the 
form W = H<frH^, where H is a Gaussian (p x n) matrix with uncorrected entries and $ 
is a (n x n) positive definite matrix that represents the covariance matrix of the channel. The 
eigenvalues distribution has been studied for the two possible cases n > p and p > n in [6], 
[7], assuming a covariance matrix $ with distinct eigenvalues (i.e., unit multiplicity). We here 
generalize the results to matrices $ with arbitrary eigenvalue multiplicities. 

Let us first recall the distributions for the case of covariance matrix with distinct eigenvalues. 



Using this result in Lemma 2 and (10) with 
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1 ) Correlation on the shortest side - distinct eigenvalues: The case p > n has been analyzed 
in [6], where it is shown that the joint p.d.f. of the (real) ordered eigenvalues Ai > A2 > • • • > A n 
of W is 

/a(xi,.,x») = — - n ;= l/f ' |V(x)| |G(x,/x)in^- n (36) 
where ^ are the n distinct eigenvalues of $ \ V(x) is the (n x n) Vandermonde matrix with 



elements Vij = x l ~ and where G(x, fx) is a (n x n) matrix with elements g it j = e~^ Xj . 

2) Correlation on the largest side - distinct eigenvalues: We here briefly derive the joint p.d.f. 
for the eigenvalues of W when <fr has all distinct eigenvalues and n > p, based on the results in 
Section III. Note that this case has been analyzed also in [7] by following a different approach. 

First we recall that, given a (p x n) random matrix H with n > p and p.d.f. 

^-png-trHHt (37) 

the p.d.f. of the (p x p) quadratic form 

W = H$H f , (38) 

where the (n x n) matrix $ is positive definite, is given by [42], [43] 

|w|"-p 

m = iH^*( r '- W )' <39) 

Then, the joint p.d.f. of the (real) ordered eigenvalues Ai > A2 > . • . > A p of W is given by 



using the results in [22, eq. (93)] as 



f x (x 1} ...,x p ) = K 1 |$r p 0^0 -W) |W|"- p Y[(xi- x 3 f , (40) 

i<j 

where 

K x = 7~~^~ • (41) 
T p {n)T p (p) 

Note that in (40) the two matrices <fr 1 and W are of dimensions (nxn) and (pxp), respectively. 
So, in (40) we evaluate F B) where B = — W © • I p is obtained by adding n — p zero 
eigenvalues to — W [7]. 

Differently from the previous literature, we can now directly use Corollary 1 and get im- 
mediately the joint p.d.f. of the ordered eigenvalues of the (p x p) matrix W when n > p 
as: 

(-l)p(n-p) TT n ,.P 

a(xi, ...,x p ) = \> n r ^ |v(x)| |G(x ' ^1 (42) 
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where fii are the eigenvalues of $ 1 , all of multiplicity one here, V(x) is the (pxp) Vandermonde 
matrix, and the (m x m) matrix G(x, fi) has elements as follows 



9i,j 



that is, the matrix G(x, fi) is 



G(x,/i) 



n—i 

Hi 



j = p + l,...,ra 



e -M*l . . . e -MiZp ^-P- 1 ^~ P ~ 2 • • • ft l 



—u,ox„ n—p—1 n—p—2 

e ^ fi 2 fi 2 



„—fj, n x p ..n—p—l ..n—p—2 



A*2 1 



fin 1 



g(*,/il) 

g(x,// 2 ) 

g(x, fin) 



(43) 



(44) 



Generalization to covariance matrix with arbitrary eigenvalues: Note that (36) and (42) 
are only valid for covariance matrices with all distinct eigenvalues (multiplicity one). So, we 
must now generalize these expressions to the case of interest, i.e., eigenvalues fii with arbitrary 
multiplicities. This step is possible by using Lemma 2. 

In fact, we note that in both (36) and (42) we have a ratio of the form 

|G(x,/x)| 



(45) 



By using Lemma 2, for each eigenvalue with multiplicity m 8 we must replace the rows of 
G(x, fi) with their successive derivatives with respect to the eigenvalue, and divide by T^^mi), 
obtaining 



|G(x sA i) 



I\i<j iv* - t*j) Hi T (md m i) Ui<j (m - m) 



niimj 



det 



r (mi-l)/ x 



g (1 W(i)) 
g ( °W(i)) 



g(^- 1 )(x,/i (L) ] 



g (0) (x,/i(L)) 



(46) 
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where the row vector g^(x, //) is the I— derivative of the row g(x, p) in (36) or (44). The j— 
element of (x, /i) is so derived to be 

(/), . I K-xjfe-^ j = l,...,p 

[ n - 3\i H n 3 3 = V + 1, • • • , n. 



The relation between the row index, i, and the derivative order, /, can be established by 
introducing the function indicating the eigenvalue fj,^ E {/^(i), • • • ^(l)} to be used in row 
i of the matrix in the RHS of (46). It is easy to verify that e» is the unique integer such that 

mi + . . . + m e ^_i < % < mi + . . . + m e .. 

Then, the derivative order for the row i is I — di, where 

di = ^2m k - i. 

k=l 

Thus, the generic element of the matrix in the RHS of (46) is gf l \jj,^). 
Combining (36), (42) and (46) we have Lemma 6. 



Appendix II 

An identity on multiple integrals involving determinants 

Theorem 2: Given an arbitrary p x p matrix <fr (x) with if— elements $i(xj), an arbitrary 
n x n matrix * (x), n > p, with elements 

^i(Xj) j = l,...,p 
*»j j =p + l,...,n 

and two arbitrary functions £(•) and £(•) the following identity holds: 



ford 



($(x)H*( x )iile(g^fe)dx 



m=l i=l 



£det({c< 

k=l ^ 



jj'=l...,n 

where the multiple integral is over the domain V ord = {b > x\ > x 2 > ■ ■ ■ > x p > a} , 

j b a ^{x)^! 3 {x) t(x)U kJ dx j = 1, . . . ,p 

j=p + l,...,n 



(48) 
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and the function Ukj(x) is denned by 



x 



if k = j 



(49) 



1 



if k^j. 



Proof: As this theorem is an extension of [6, Theorem 3], it is sufficient for the proof to 
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Fig. 2. Ergodic mutual information for single-user MIMO systems as a function of SNR over Rayleigh uncorrelated fading 
with m = 6, riR = 3. Half of the antennas with power (normalized) 1 + A, the others with 1 — A, i.e., with transmitted power 
levels (normalized) equal to {1 + A, 1 + A, 1 + A, 1 - A, 1 - A, 1 - A}. 
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Fig. 3. Bound on the network capacity for MIMO relay networks. Source with 4 antennas, 5 relays with 2 antennas each, 
power levels per relay proportional to {1, 2, 5, 10, 20}. 
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Fig. 4. Ergodic mutual information for MIMO-(6,6) as a function of SIR in the presence of one MIMO cochannel interferer 
with Nil — 1, 2, 4, 6, 10. The SNR is set to 10 dB. The Gaussian approximation of the interference is also shown. Diamond: 
capacity of a single-user MIMO-(6, 6). Circles: capacity of a single-user MIMO-(6, 6 — iVri) (only for jVri = 1, 2, 4). 
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SIR (dB) 



Fig. 5. Ergodic mutual information as a function of the signal-to-total interference ratio. MIMO system with TVr = 6 receiving 
antenna, SNR = 10 dB. The Gaussian approximation of the interference is also shown. Scenario with one and two interferers, 
each with the same number of transmitting antennas as the desired user. Cases of 3, 4, 5 and 6 transmitting antennas. Circles: 
capacity of single-user MIMO-(A r To, Nr). 



DRAFT 



30 



SUBM. TO IEEE TRANS. ON INF. TH. 



Marco Chiani (M'94-SM'02) was born in Rimini, Italy, in April 1964. He received the Dr. Ing. degree 
(magna cum laude) in Electronic Engineering and the Ph.D. degree in Electronic and Computer Science 
from the University of Bologna in 1989 and 1993, respectively. Dr. Chiani is a Full Professor at the II 
Engineering Faculty, University of Bologna, Italy, where he is the Chair in Telecommunication. During the 
summer of 2001 he was a Visiting Scientist at AT&T Research Laboratories in Middletown, NJ. He is a 
frequent visitor at the Massachusetts Institute of Technology (MIT), where he presently holds a Research 
Affiliate appointment. 

Dr. Chiani's research interests include wireless communication systems, MIMO systems, wireless multimedia, low density 
parity check codes (LDPCC) and UWB. He is leading the research unit of University of Bologna on cognitive radio and 
UWB (European project EUWB), on Joint Source and Channel Coding for wireless video (European projects Phoenix-FP6 and 
Optimix-FP7), and is a consultant to the European Space Agency (ESA-ESOC) for the design and evaluation of error correcting 
codes based on LDPCC for space CCSDS applications. 

Dr. Chiani has chaired, organized sessions and served on the Technical Program Committees at several IEEE International 
Conferences. In January 2006 he received the ICNEWS award "For Fundamental Contributions to the Theory and Practice of 
Wireless Communications". He was the recipient of the 2008 IEEE ComSoc Radio Communications Committee Outstanding 
Service Award. 

He is the past chair (2002-2004) of the Radio Communications Committee of the IEEE Communication Society and past 
Editor of Wireless Communication (2000-2007) for the IEEE TRANSACTIONS ON COMMUNICATIONS. 



PLACE 
PHOTO 
HERE 



DRAFT 



CHIANI, WIN, SHIN: MIMO NETWORKS. 



31 



Moe Z. Win (S'85-M'87-SM'97-F'04) received both the Ph.D. in Electrical Engineering and M.S. in 
Applied Mathematics as a Presidential Fellow at the University of Southern California (USC) in 1998. He 
received an M.S. in Electrical Engineering from USC in 1989, and a B.S. (magna cum laude) in Electrical 
Engineering from Texas A&M University in 1987. 

Dr. Win is an Associate Professor at the Massachusetts Institute of Technology (MIT). Prior to joining 
MIT, he was at AT&T Research Laboratories for five years and at the Jet Propulsion Laboratory for seven 
years. His research encompasses developing fundamental theories, designing algorithms, and conducting experimentation for a 
broad range of real-world problems. His current research topics include location-aware networks, time-varying channels, multiple 
antenna systems, ultra-wide bandwidth systems, optical transmission systems, and space communications systems. 

Professor Win is an IEEE Distinguished Lecturer and elected Fellow of the IEEE, cited for "contributions to wideband 
wireless transmission." He was honored with the IEEE Eric E. Sumner Award (2006), an IEEE Technical Field Award for 
"pioneering contributions to ultra-wide band communications science and technology." Together with students and colleagues, 
his papers have received several awards including the IEEE Communications Society's Guglielmo Marconi Best Paper Award 
(2008) and the IEEE Antennas and Propagation Society's Sergei A. Schelkunoff Transactions Prize Paper Award (2003). His 
other recognitions include the Laurea Honoris Causa from the University of Ferrara, Italy (2008), the Technical Recognition 
Award of the IEEE ComSoc Radio Communications Committee (2008), Wireless Educator of the Year Award (2007), the 
Fulbright Foundation Senior Scholar Lecturing and Research Fellowship (2004), the U.S. Presidential Early Career Award for 
Scientists and Engineers (2004), the AIAA Young Aerospace Engineer of the Year (2004), and the Office of Naval Research 
Young Investigator Award (2003). 

Professor Win has been actively involved in organizing and chairing a number of international conferences. He served as the 
Technical Program Chair for the IEEE Wireless Communications and Networking Conference in 2009, the IEEE Conference on 
Ultra Wideband in 2006, the IEEE Communication Theory Symposia of ICC-2004 and Globecom-2000, and the IEEE Conference 
on Ultra Wideband Systems and Technologies in 2002; Technical Program Vice-Chair for the IEEE International Conference on 
Communications in 2002; and the Tutorial Chair for ICC-2009 and the IEEE Semiannual International Vehicular Technology 
Conference in Fall 2001. He was the chair (2004-2006) and secretary (2002-2004) for the Radio Communications Committee of 
the IEEE Communications Society. Dr. Win is currently an Editor for IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS. 
He served as Area Editor for Modulation and Signal Design (2003-2006), Editor for Wideband Wireless and Diversity (2003- 
2006), and Editor for Equalization and Diversity (1998-2003), all for the IEEE TRANSACTIONS ON COMMUNICATIONS. He 
was Guest-Editor for the PROCEEDINGS OF THE IEEE (Special Issue on UWB Technology & Emerging Applications) in 2009 
and IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS (Special Issue on Ultra -Wideband Radio in Multiaccess 
Wireless Communications) in 2002. 



PLACE 
PHOTO 
HERE 



DRAFT 



32 



SUBM. TO IEEE TRANS. ON INF. TH. 



Hyundong Shin (S'Ol-M'04) received the B.S. degree in Electronics Engineering from Kyung Hee 
University, Korea, in 1999, and the M.S. and Ph.D. degrees in Electrical Engineering from Seoul National 
University, Seoul, Korea, in 2001 and 2004, respectively. 

From September 2004 to February 2006, Dr. Shin was a Postdoctoral Associate at the Laboratory 
for Information and Decision Systems (LIDS), Massachusetts Institute of Technology (MIT), Cambridge, 
MA, USA. In March 2006, he joined the faculty of the School of Electronics and Information, Kyung 
Hee University, Korea, where he is now an Assistant Professor at the Department of Electronics and Radio Engineering. His 
research interests include wireless communications, information and coding theory, cooperative/collaborative communications, 
and multiple-antenna wireless communication systems and networks. 

Professor Shin served as a member of the Technical Program Committee in the IEEE International Conference on Communica- 
tions (2006, 2009), the IEEE International Conference on Ultra Wideband (2006), the IEEE Global Communications Conference 
(2009, 2010), the IEEE Vehicular Technology Conference (2009 Fall, 2010 Spring), and the IEEE International Symposium 
on Personal, Indoor and Mobile Communications (2009). He served as a Technical Program co-chair for the IEEE Wireless 
Communications and Networking Conference PHY Track (2009). Dr. Shin is currently an Editor for IEEE TRANSACTIONS 
on Wireless Communications. He was a Guest Editor for the 2008 EURASIP Journal on Advances in Signal 
PROCESSING (Special Issue on Wireless Cooperative Networks). 

Professor Shin received the IEEE Communications Society's Guglielmo Marconi Prize Paper Award (2008) and the IEEE 
Vehicular Technology Conference Best Paper Award (2008 Spring). 



PLACE 
PHOTO 
HERE 



DRAFT 



